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The flow and deposition of polydisperse granular materials is simulated through the Magnetic 
Diffusion Limited Aggregation (MDLA) model. The random walk undergone by an entity in the 
MDLA model is modified such that the trajectories are ballistic in nature, leading to a magnetically 
controlled ballistic deposition (MBD) model. This allows to obtain important ingredients about 
a difficult problem that of the nonequilibrium segregation of polydisperse sandpiles and heteroge- 
neous adsorption of a binary distribution of particles which can interact with each other and with 
an external field. Our detailed results from many simulations of MBD clusters on a two dimensional 
triangular lattice above a flat surface in a vertical finite size box for binary systems indicates intrigu- 
ing variations of the density, "magnetization", types of clusters, and fractal dimensions. We derive 
the field and grain interaction dependent susceptibility and compressibility. We deduce a completely 
new phase diagram for binary granular piles and discuss its complexity inherent to different grain 
competition and cluster growth probabilities. 



I. INTRODUCTION 

Understanding the flow and static structures of granu- 
lar matter is becoming increasingly relevant. Many raised 
questions are tackled along various lines of approach[l| . 
It has already been claimed that the simplifications found 
in basic models need to be improved in order to explain 
features of such complex materials. Granular pile spread- 
ing processes driven by cooperative non-linear evolution 
rules lead to developed patterns which often reach a high 
level of complexity. It is of present interest to examine 
whether growth models can be used for describing gran- 
ular structures and related material properties. Further- 
more the non-linear processes at work in granular flows 
and depositions hint toward simulation approaches^- 
Cooperative effects in ballistic deposition of hard disks 
have been recently mentioned 3]. 

Kinetic growth models (KGM) have received much 
attention, |j| like the Eden model[5j and the diffusion lim- 
ited aggregation model@ (DLA). They have served to 
describe nonequilibrium phenomena like film or crystal 
growth0,Q, epidemics^, material fractures [Tof. ... In 
all cases, such models are mainly concerned by the tran- 
sition from dense branching to dendritic morphology. 

One important physical constraint has to be considered 
in describing granular materials: the materials are not 
made of symmetrical (spherical or cubic) entities. The 
surface of grains is usually rough, thus leading to spe- 
cific angles of repose 11]. Also the grain anisotropy leads 
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to phenomena like iams[T^|. in flow, and arches m 
static structures. It seems therefore necessary to have at 
least one degree of freedom in order to describe grains; 
we are even aware that only one degree is a very strong 
approximation. This degree of freedom should be cou- 
pled to some field, just like a spin to a magnetic field. 
Whence one can imagine that grains are identical enti- 
ties except for one degree of freedom, call it a spin though 
it can be any physical feature of particular interest, like 
the grain roughness or shape feature. Clearly a spin al- 
lows for referring to a direction or a rotation process; if 
this is admitted, to take such a degree of freedom into 
account in describing granular piles should basically im- 
prove the granular state overall description. (Generaliza- 
tions are immediately imagined by anyone familiar with 
spin models and statistical mechanics; one can later on 
imagine many component vector models, including Potts- 
like models |lj|). In fact a constrained Ising spin chain 
has been recently considered and studied as a toy model 
for granular compaction|irij . The exchange energy J de- 
scribing the " spin-spin interaction" is analogous in granu- 
lar matter to the contact energy due to surface roughness 
between grains. A similar interpretation of J for flows 
can be found in Pandey et al.Q. The external "mag- 
netic" field in such a case can be e.g. a wind field, the 
sign depending on e.g. change in pressure due to grain 
drag. We thus combine topology and mass (or weight) in 
order to describe granular materials in a simple way. 

The above ideas remind us that a similar set of 
considerations has been found in the magnetic Eden 
model^j and in the Ma gne tic Diffusion Limited Ag- 
gregation (MDLA) model[l7j when attempting to de- 
scribe crystal growth in a magnetic field, when there is 
a competition between entities. Aggregation can pro- 
ceed under short range or long range dipolar interactions 
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in fact ll9j- However the studies pertained to the 
growth of clusters starting from a point seed. In a re- 
cent set of investigations on the magnetic Eden model, 
Albano et al.^HjllilEi have pointed out the interest 
of such models for examining deposition and film growth, 
thus starting from a substrate, - sometimes with rather 
complex realistic rules. The studies are also related to 
nonequilibrium wetting questions. |24| and other deposi- 
tion problems, plf 

The same type of studies can be done with the MDLA, 
i.e. examining the growth of clusters from a substrate. 
It is clear that a substrate having finite (or not) size 
destroys the spatial isotropy or more generally spatial 
symmetry. In simple words, there is a top and bottom, 
if the substrate is horizontal. It is obviously of interest 
for granular materials in a gravitational field to consider 
what happens only in the half space above the (finite 
size of course) substrate. This reduction is however of 
fundamental and practical interest because one can also 
consider that the system is in a vertical box with walls, as 
in the Albano et al. film growth geometry [21], with 
binary competing entities. In DLA, the diffusing particle 
follows a random walk.Q However, it seems very hard to 
let this usual DLA rule holds here concerning the path 
of the granular entity launched far away from the seed or 
substrate. In the present considerations, it seems more 
appropriate to let the granular entity follows a ballistic 
vertical trajectory like in rain modelsfHJ rather than a 
random walk. 

In this paper, we calculate what changes result in the 
features of a classical ballistic deposition |25l I26L (BD) 
model when we add one extra degree of freedom, a " spin" , 
to the classical BD. The spin reflects e.g. the orientation 
of oblate grains or their roughness characteristic. The 
spins, as usual, interact through some exchange energy 
J which, for grains, is often mechanical or electrostatic 
in nature. The external field H is thought to be the 
image of a classical field positioning or influencing the 
flow of grains, like some wind velocity or, more gener- 
ally, pressure difference, or an electric field. This model 
is hereby called the magnetically controlled ballistic de- 
position (MBD) model falls into the category of kinetic 
growth models. In Sect. 2, we enumerate the algorithm 
rules and briefly comment upon them. In Sect. 3 we 
present results on the "density" and "magnetization", 
"susceptibility" and " compressibility" . The types of clus- 
ters and their fractal dimensions are discussed in Sect. 4. 
In Sect. 5, a brief conclusion can be found. 



II. EXPERIMENTAL PROCEDURE 

For simplicity we will thereafter call a grain, a spin. 
It will take here only two states (up or down) or two 
values (+1 or -1). The external field is supposed to be 
constant and uniform throughout the whole system. For 
obvious reasons, like higher packing considerations and 
possible geometrical frustrations, the underlying lattice 
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FIG. 1: Example of a 4x4 triangular lattice. The arrow 
indicates the fall direction. The substrate has been marked 
with crosses. 

should not have a square symmetry; we have taken a 
triangular lattice in the following simulations (Fig^J 

One of the main problems to be tackled is how to 
choose the best rule for aggregating granular falling enti- 
ties (spins) in order to get them stick together and form 
clusters. We have chosen the simplest sticking rule, usual 
in statistical physics, i.e., namely a Metropolis-like rule, 
as also in Pandey et al.Q, presented here below. The 
algorithm goes as follows: 

1. first, we choose a horizontal substrate of spins with 
a predetermined (for example antiferromagnetic- 
like or random) configuration; 

2. a falling (up or down) spin is dropped along one of 
the lattice lines (see Fig^l from a height r max + 5a, 
where r max is the largest distance between a cluster 
site and the substrate i.e. here it is the height of 
the highest column growing from the substrate on 
the lattice; at each step down the spin can flip i.e., 
change its sign, with equal probability; e.g. the 
anisotropic grain can rotate; 

3. the spin goes down until it reaches a site perimeter 
of the cluster; the local gain in the Ising energy 

<i.j> i 

is calculated before and after the spin possible im- 
pact, thus cluster growth. If the gain is negative 
the spin sticks to the cluster immediately (stick- 
ing probability =1.0) and we go back to step 
In the opposite case the spin sticks to the cluster 
with a rate exp(— A(3E) where A/3E is the local 
gain in the Ising energy. If the spin does not stick 
to the cluster it continues going down toward the 
substrate or bottom of the box. Of course if the 
site just below the spin is occupied the spin imme- 
diately stops and sticks to the cluster. When the 
spin sticks to the cluster we go back to step (J2J). 

4. After dropping a (large) number of spins the phys- 
ical quantities of interest like the magnetization, 
density, fractal dimension, etc. are computed. 
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FIG. 2: Relations between J and j3H on the f5J-f3H plane; 
32 regions result from these relations. 



It should be noticed that there is no toppling nor re- 
laxation at this time like in Manna sandpile model or 
its extensions'^ l28l l29l |. Moreover, since the number 
of nearest neighbors on a triangular lattice is equal to 
6, and due to the rule of MBD (above paragraph), we 
do not have to take into considerations the spin configu- 
rations in which the depositing spin has a neighbor just 
vertically over it or just below (because it would then 
always stick to the cluster). Therefore there are only 4 
neighboring sites where spin configurations are relevant 
for calculating the local gain in the Ising energy. On the 
other hand there are 3 kinds of site occupation : spin up, 
spin down and no spin, hence 52 configurations (exclud- 
ing the empty one - when there is no spin on the perime- 
ter). Some of these configurations are symmetrical with 
respect to rotations. Finally 23 configurations are to be 
examined having at least one spin on the perimeter. All 
these sticking configuration rates (also distinguishing the 
sign of the falling spin) are shown in table [I] The contri- 
bution to the sticking rate arising from the interaction of 
the depositing spin with the field has also to be evaluated. 

The rates of sticking to the cluster read like: 



P n „ = e~ Al3E = ^PJ+PH)^ (2) 

where Af3E is the local gain of the Ising energy, and 
n is the difference between the number of up and down 
spins: the possible values are —4, —3, ..,4, and s is the 
value (or sign) of the falling spin(— 1 or 1). Equating all 
these rates lead us to a set of 16 relations between (3J 
and (3H; in fact as in the MDLAp^]. FigEl determine 32 
regions where granular packing cluster growth processes 
differ from each other. The case (3J = and j3H = 
corresponds to the usual ballistic deposition model (stan- 
dard deposition without spin and field). y 

Notice that the model contains two order parameters, 
as in the Blume-Emery-Griffiths model[3(], one cor- 
responds to the density, the other to the magnetization. 
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FIG. 3: Dependence of density on f3H for different /3J. 
Observe that the behavior is symmetrical with respect to 
f3H = 0. 

III. NUMERICAL RESULTS AND DISCUSSION 

All results reported below are for a triangular lattice 
of horizontal size L = 100, i.e. the width of the seed sub- 
strate, and when the pile made of clusters has reached a 
500 lattice unit height. Every reported data point corre- 
sponds to an average over 1000 simulations. 

A. Density 

We define the density of a cluster as 

q number of spins in the cluster 
number of sites on the lattice ' 

in which obviously the number of lattice sites = 50000. 
FigO illustrates the behavior of the density with respect 
to the (3H parameter. This figure convinces us that the 
results are symmetrical with respect to (3H = 0. There- 
fore in the following subsections we will often present 
results for f3H > only. It is observed that the density 
presents a sharp minimum when j3H = ±(3J. The gran- 
ular pile is rather loosely packed since the density varies 
between 0.37 and 0.45. This low value with respect to 
experimental findings arises from the fact that we have 
not included relaxation processes in this investigation. 

On figB] one illustrates in a 3D way the influence of 
the interaction part of the hamiltonian, i.e. f3J, on the 
density. This sort of diagram allows us to emphasize 
that the density is almost the same everywhere, as men- 
tioned above, but there are several minima: the main 
density variations occur along lines bordering plateau re- 
gions, lines which correspond to the equal probability 
lines mentioned in the previous section. One can ob- 
serve a set of trenches near these borders between differ- 
ent growth (or packing) regions. In these trenches, the 
density is markedly lower than in the immediate neigh- 
borhood. Indeed such trenches correspond to the highest 
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TABLE I: Rate of sticking and configurations in the MDM model on the triangular lattice. First (second) column shows the 
rate of sticking an up (down) spin to the preexisting spin configuration. 



possible rate of sticking, thus to a condition for loose 
packing (see also below). 
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FIG. 5: Dependence of the density on pV for different f3H . 



FIG. 4: Dependence of the MBD pile density on [5H and f3J . 
Observe trenches at borders between plateau regions. 

From Fig03 the value of the density is seen to be 
slightly different for f3J > and for f3J < 0, but re- 
mains qualitatively the same. These differences indicate 
that a little higher density is obtained for ferromagnetic 
systems f3J > than for antiferromagnetic ones, in par- 
ticular when the external field is different from 0. Some- 
thing similar had been found in studies on the MDLA[I3]. 
Further discussion on this point is postponed for after ex- 



amination of the clusters in Sect. 4. 

Another illustration of the density dependence is ex- 
hibited on FigEJas a projection on the (/3J, (5H) plane. 
The main trenches are observed, i.e. only 10, out of all 
32 possible ones: they are located at f}H — 0, (3H = 
±/3J, j3H — ±2/3,7. One can distinguish a trench for 
positive and negative values of f3J and (3H . The last 
one is not clearly visible on FigEl but on Fig[5] one 
can observe a small hollow near the trench j3H = —(3J 
for (3J < 0. Hollows are positioned along the trench 
f3 = — 2(3J. Such structures (hollows) are artefacts due 
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FIG. 6: Dependence of the density on (3H and (3J . Different 
gray colors correspond to different density levels with "color" 
scale indicated on the right. 

to resolution of the simulation. 

The (3H = ±2(33 trench is very clearly seen on FigJSJ 
for (3H = 5 a small minimum is observed for (33 = ±5, a 
symptom of the existence of a trench. 

To explain the behavior of the density due to (33 and 
(3H parameters let us compute the maximum probability 
of cluster sticking for the configuration having n spins in 
a site neighborhood. The results for each configuration, 
with respect to the number of spins in the perimeter arc 
presented on Fig[7| The discrete color changes is real but 
the zig zag blurring is an artefact due to the numerical 
sampling with (33 and (3H . Next, from this figure it can 
be concluded that the biggest drop in density occurs for 
(3H = ±/3 J indeed (FigEK)- The reason stems from the 
configurations which correspond to a case with a small 
amount of spins i.e. with 1 or 3 spins, thus to a possibly 
high relative change in Ising energy. This observation 
further explains the existing density trenches observed 
on FigH (or Fig© for (3H = ±/3 3. 

A very similar conclusion pertaining to the (3H = 
and (3H = ±2/3 3 lines can be drawn for the appropriate 
trenches in Fig0| (or Fig|SJ); the drop is thereby smaller 
because these configurations contain 2 and 4 spins, - 
starting from a 4 spin configuration does not obviously 
reduce the density after one extra spin sticks to the clus- 
ter. 

At each step, the available "volume of interest" on the 
perimeter corresponds to 7 sites. For a configuration with 
the largest possible number of "useful" spins for stick- 
ing in the perimeter we have a final 5 spin configuration 
(4 + the dropping spin), whence the local density is at 
most | w 0.71. However this occurs at most 6 times 
(out of 23), see Table I leading to a rough estimate of 
2(6/23)(5/7)=0.37 for the density. 

On the other hand, configurations with 3 and 4 spins 
on the perimeter, lead to an increase in the density after 
spin sticking roughly equal to 0.46 as observed. 

The behavior of the density in the vicinity of the zero 



field and zero interaction case is of interest with respect 
to observe deviations from the classical BD model. On 
FiglHl it is seen that when the field is different from 
the density is almost constant over some (3 3 interval, the 
width of which depending on the field value, i.e. the 
higher the field, the wider the interval. Observe that the 
density value in the interval does not seem to be varying 
linearly with the field. Notice that the density depen- 
dence observed on Fig0and its characteristic structures, 
i.e. the trenches are again observed for moderately high 
fields (ca. 0.8). 

Finally, FigEl exhibits the behavior of the density for 
f3H = (33 i.e., along the deepest trench (FigQJor Fig|HJ). 
It should be emphasized that the (3 3 < trench level has 
a lower density than the (3 3 > one due to the imbal- 
ance in sticking probabilities for preferred ferromagnetic 
or antiferromagnetic-like configurations. Worth mention- 
ing is that the second deepest trench, i.e. (3H = —j3H has 
exactly the same dependence like that shown on Fig|5| 

B. Magnetization 

In this subsection we present results concerning the 
magnetization of the cluster of spins. The magnetization 
is defined as: 



where n + and rt_ are the number of up and down spins 
respectively. This quantity can be considered as a mea- 
sure of the difference in grain orientations in the packing. 

As seen on Fig^| the magnetization dependence on 
(3H and (33 presents a sort of terraces and is slightly 
undulating, at the borders of specific regions previously 
emphasized in the density dependence discussion. Inter- 
estingly all expected dependences are better visible on 
the magnetization than on the density pictures. In par- 
ticular, see Figllll 

Figll2l and Fig^| present the behavior of the magne- 
tization for different (3J and (3H values. The maximum 
magnetization occurs for (3 J = and no field. For a finite 
field the maximum is rather broad. Terrace structures 
are seen like in Fig^Jfor the density. There are about 6 
terraces with different values of the magnetization. Each 
terrace occurs in ranges like those of the regions previ- 
ously mentioned. It is emphasized that the magnetiza- 
tion terrace levels differs in the antiferromagnetic and 
the ferromagnetic coupling regions, i.e. for (33 < and 
(33 > 0, - the level height depending on the field sign. 

This is well stressed through figure fFig ll3l) which ex- 
hibits the dependence of the magnetization on (3H for 
different (33. Notice that for the case without spin-spin 
interaction, one obtains a kind of saturation at a value 
ca. ±0.36. Observe the surprising form of the (M, H) 
curve, as for classical soft magnets. For a finite field, the 
magnetization does not saturate like in the case of zero 
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FIG. 7: Theoretical maximum probability of sticking a spin to a cluster for a cluster configuration with: (a) 1 spin, (b) 2 spins, 
(c) 3 spins, (d) 4 spins. 
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FIG. 8: Dependence of the density on f3.J for different (3H. 



FIG. 9: Dependency of the density on f3J for (5H — /3J i.e., 
the deepest trench. 



field, for the values that we have investigated, but M 
saturates creating stair-like structures. 



C. Magnetic susceptibility 

The magnetic susceptibility of the clusters can be ob- 
tained by numerically differentiating the magnetization 
M over the field H , at fixed J or H , i.e. 



XJ = 



dM 



XH = 



dM 
~dJ 



(5) 



H 



The dependence is presented on FigEl The highest 
susceptibility occurs for regions where (3H — (3J = 
and the main trenches. Other regions have a rather rela- 
tively small susceptibility. The results can be understood 
through the role of the interaction between spins which 
as usual induces a drop in the magnetic susceptibility 



FIG. 10: Dependence of the magnetization on f3J and f3H . 



FIG. 13: Dependence of the magnetization on f3H for differ- 
ent PJ values. 





- - - 0.00 



FIG. 11: Dependence of the magnetization on /3J for different 
PH, with a gray color scale as indicated. 
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FIG. 12: Dependence of the magnetization on pj for different 
pH. 



FIG. 14: Dependence of the magnetic susceptibility on pj 
and PH. 



of the materials. On the other hand, the variation of 
the orientation difference (M) of the grains in a packing 
with respect to an external (wind or) field should be an 
interesting experimental test. 



D. Compressibility 

Similarly the compressibility (Figll5|) at fixed H or J, 

is 



1 dG 
G~dJ 



H 



1 dG 
GdH 



(6) 



where G is the density. The displayed data is rather 
blurred because of the limited amount of data for nu- 
merical differentiation near the trenches and in the re- 
gion of the standard BD, in particular. Nevertheless this 
figure indicates some mild variation due to internal com- 
petition and external conditions. Let us also recall that 
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FIG. 15: Dependence of the compressibility on /3J and /3H. 

while the susceptibility is singular in spin glass models 
of compaction [32j at the critical percolation value, the 
compressibility seems to remain finite. 

0,113 Our model 

indicates the same, as experimentally or numerically ob- 
served. Further experimental and numerical considera- 
tions should be given to this point. 

IV. PILE STRUCTURES 

In this section we present some examples of typical 
clusters created by the MBD in specific regions, as ob- 
served and discussed here above. The size of the lacunes 
in each cluster allows some emphasis and contrasting. 



the internal competition leads to different cooperative 
phenomena in cluster packing. [3j 



B. Fractal dimension 

For further relating the model and our investigations 
to granular piles, it is of interest to check the fractal 
dimension yj, |35j of the piles in the different parame- 
ter regions. The (box counting) technique |34l l35| con- 
sists in covering, without overlapping, the whole cluster 
by squares of the same size, and computing the number 
of squares which have at least one spin up (or down), 
for different square sizes. We have distinguished be- 
tween the fractal dimension of the cluster of up and 
down spins. When computing the fractal dimension of 
the whole cluster it was simply checked whether there 
was at least one spin and its sign in the relevant square. 
From the best linear fit to the data, i.e. — log(square size) 
vs. log(number of squares with a spin), the fractal di- 
mension is obtained through the slope. The results are 
reported in Figs 1 171 181 

Let us discuss the results corresponding to differ- 
ent regions, either not distinguishing over the spin sign 
fFig ll7f) . or on the contrary considering down or up cases 
fFiglTHh or Fig^Js, respectively). The fractal dimension 
in every case ranges from 1.91 to 1.95, therefore is about 
equal to 2, taking into consideration the error bars. This 
value is similar to what is found in classical BDM0,Ei 
and in the rain model |25j. Surprisingly the lowest fractal 
dimensions are found in the above for the trenches and 
in particular for the classical BDM. 



A. Typical clusters 

In Fig llfil 9 clusters from 9 different growth regions 
are shown . The central cluster has the smallest density: 
it corresponds to the case when there is no interaction 
between spins and no field, i.e. it is the standard BD 
model j2g. 

On the other hand the clusters for (3J = and 
f3H = ±5 have the highest density; they illustrate regions 
where the density is saturating (above 0.44, see Fig|3J|. 
It should be mentioned that the results are symmetrical 
with respect to the field sign. 

For a finite J interaction between spins, a quite dif- 
ferent behavior of the clusters is observed. The clusters 
corresponding to the trenches (j3J = ±5 and (3H = ±5) 
have a smaller density (ca.0.41) than those outside the 
trenches (see for example clusters with (3 J = and 
f3H = ±5). Also observe, in Fig lltil field-free grown clus- 
ters and differences between (3 J < and (3 J > clusters. 
When (3, J > 0, the spins show a tendency toward similar 
sign spin " domains" . In the antifcrromagnetic-like re- 
gion, i.e. where f3J < 0, adjacent spins have more often 
opposite directions, - a cluster with (3 J < and (3H = 
is a typical example. This allows us to emphasize that 



V. CONCLUSION 

We have presented a nonequilibrium ballistic deposi- 
tion model with one degree of freedom per entity, degree 
which can be coupled to an external field. We have ex- 
amined the cluster properties emphasizing the existence 
of two order parameters, since two characteristic field (J 
and H) are intrinsic to the model. This model can serve 
to describe in a first approximation the deposition of a 
distribution of grains, distribution characterized by one 
intrinsic parameter which can be coupled to an external 
field. The degree of freedom can be either the anisotropy 
factor or the surface roughness or an electrostatic imbal- 
ance of a grain, - the corresponding field being immedi- 
ately thought of. 

For the sake of such an extension of usual deposition 
models, the degree of freedom has been called "spin". 
We have simulated the nonequilibrium deposition in a 
finite size 2D vertical box, admitting that grains flow 
down along linear trajectories on a triangular lattice. The 
" quenching" of the degree of freedom on the cluster leads 
to branching or compactness and moreover to combined 
geometric and physical regions at specific "field" and 
"spin-spin interaction" values. This was seen through the 
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FIG. 16: Typical clusters. The "colors" correspond to the spin's sign. 




FIG. 17: The fractal dimension for the whole cluster, i.e. not 
recognizing the spin signs. 

calculation of the " density" and the so called " magnetiza- 
tion" . Different cluster regimes were expected according 



to the (spin sticking) packing rule. 

We have investigated a box geometry and a triangular 
underlying lattice. The ranges in density and magnetiza- 
tion are limited, but features exist resulting from compet- 
itive nonequilibrium growth/deposition processes. Min- 
ima in density occur along specific sticking probability 
lines. Slight differences exist whether the "spin-spin in- 
teraction energy" is positive or negative. The fractal di- 
mension of clusters whatever the type of grains and the 
parameter sign or values is however rather trivial and 
equal to 2. This differs markedly from what was found in 
the MDLA model starting from a central seed in which 
both in the ferromagnetic interaction regions, and the 
AF regions, the cluster morphology was dendritic with 
an important thickening of the branches and the fractal 
dimension ranging from 1.68 to 1.99. Instead of a criti- 
cal a value /3J C a set of values depending on the external 
field divide the parameter plane in regions, with trenches 
and plateaus. Thus, a spreading phenomenon is avoided 
around such f3J c (H). It is of interest to further exam- 
ine whether this has interesting consequences in granular 



10 





FIG. 18: The fractal dimension for the granular cluster when we recognize the spin signs. 



deposition situations. |3fij 

The connection between this model and granular mat- 
ter systems suggests some experimental work. In the 
MBD the spin can be interpreted as a rotation or defin- 
ing a direction process. The coupling constant can be 
mapped to a mechanical friction energy, the magnetic 
field to gravity or wind pressure. It is known that pack- 
ing is impaired by static electrical charges. Such effects 



may be considered in the above framework. 
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